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Abstract. - We present an exact solution describing equilibrium properties of the catalytically- 
activated A + A — > reaction taking place on a one- dimensional lattice, where some of the 
sites possess special "catalytic" properties. The A particles undergo continuous exchanges with 
the vapor phase; two neighboring adsorbed As react when at least one of them resides on a 
catalytic site (CS). We consider three situations for the CS distribution: regular, annealed 
random and quenched random. For all three CS distribution types, we derive exact results 
for the disorder-averaged pressure and present exact asymptotic expressions for the particles' 
mean density. 



Introduction. - Catalytically-activated reactions (CARs), i.e. reactions involving par- 
ticles which may recombine only in the presence of some third substance - a catalyst - are 
widespread in nature pQ. The properties of such reactions have attracted considerable interest 
recently, following an early analysis of a particular reaction scheme, namely the CO-oxidation 
in the presence of metal surfaces with catalytic properties It has been realized that this 
reaction exhibits an essentially different behavior when compared to the predictions of the clas- 
sical, formal kinetic scheme pQ and that under certain conditions such collective phenomena 
as phase transitions or the formation of bifurcation patterns may occur [SJIIIIE]. 

A common assumption of these studies 12131111 is that the catalyst is modelled as an ideal 
surface with homogeneous catalytic properties. On the other hand, in realistic systems, the 
catalyst is often not a well-defined object, but rather consists of mobile or localized catalytic 
sites (CSs) or islands, whose spatial distribution is complex £Q. Metallic catalysts, for in- 
stance, are often disordered compact aggregates, the building blocks of which are imperfect 
crystallites with broken facets, kinks and steps. In porous materials with convoluted surfaces, 
such as e.g., silica, alumina or carbons, the effective catalyst occupies only a portion of the 
total surface; in solution, the catalyst can consist of active groups attached to polymer chains. 

© EDP Sciences 



2 



EUROPHYSICS LETTERS 



• • VAPOR • ^ 

(a) 

i % i T T x % T x X X 1 

1 N N+1 

Fig. 1 - ID lattice of N adsorption sites in contact with the vapor phase. The filled circles denote 
A particles. The crosses denote the adsorption sites with catalytic properties (CS). (a) denotes a 
"forbidden" particle configuration. 



Such complex morphologies render the theoretical analysis difficult and, as yet, only em- 
pirical approaches have been proposed (see, e.g. Ref. |TJ). Consequently, analytical solutions 
even of idealized or simplified models are highly desirable, since such studies may provide an 
understanding of the effects of heterogeneities on the properties of CARs. 

In this Letter we study the equilibrium properties of a heterogeneous, catalytically-activated 
A + A — > reaction in a simple, one-dimensional (ID) model. The catalyst is modelled here 
as an array of special catalytic sites (CSs). In regard to the CS distribution, we focus on 
three situations: regular, annealed random and quenched random. For all three cases we 
derive exact results, taking into account equilibrium fluctuations, for the disorder-averaged 
pressure and present exact asymptotic expressions for the particles' mean density. We show 
that despite the apparent simplicity of the model, one obtains highly non-trivial behaviors. 

Model. - Consider a ID regular lattice containing TV adsorption sites (Fig.l), which is 
brought in contact with a reservoir (vapor phase) of non-interacting A particles with hard- 
cores, maintained at a constant chemical potential /i; the activity is hence z — exp(/3/x). 

The A particles from the vapor phase can adsorb onto vacant adsorption sites and desorb 
back to the reservoir. The occupation of the "i"-th adsorption site is described by the variable 
m such that ni = 1 if the "i"-th site is occupied, and rii = otherwise. We assume that some 
of the adsorption sites possess "catalytic" properties (denoted by crosses in Fig.l) in the sense 
that they induce an immediate reaction A + A — ► between neighboring A particles; that 
is, if at least one of two neighboring adsorbed A particles resides on a catalytic site, these 
two particles instantaneously react and leave the chain. In case when both neighboring sites 
are occupied, the A particle landing on the CS reacts with either of its neighbors with equal 
probability. The positions of the catalytic sites are specified by the variable Q, such that Q = 1 
if the "i"-th site is catalytic, and Q — otherwise. For convenience, we add two boundary 
sites, i.e. i = and i = N+1, and stipulate that these sites are always non-catalytic, i.e. 
Co = Cn+i = 0, and always unoccupied, no = rtw+i = 0. 

We focus here on the equilibrium properties of the model under study. For a given CS 
distribution, the partition function Zn{Q of the- adsorbate, being in equilibrium with the 
reservoir and constrained locally by the condition that no two A particles can occupy simul- 
taneously two neighboring sites if at least one of them is a CS, can be written as: 

N 

Zn(C) = II z "' ( x - (im n i-i) (l - dniUi+i) . (1) 

{n ^ }^=l 
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In this way, any two neighboring sites i and i — 1 are coupled by the factor (1 — mrii-i) 
when at least one of these sites is catalytic. These coupling factors are depicted in Fig.l as 
arcs connecting neighboring sites; note that the configurations {rii} in which the variables 
and Ui-i of two connected sites assume both the value 1 are excluded. Note also that the 
partition function in eq.([Q) involves effectively three-site interactions. We introduce now the 
notion of a " iC-cluster" , as being a set of K sites connected to each other by arcs. Note that 
the boundary between adjacent clusters is given by a pair of two neighboring, non-catalytic 
sites, i.e. when two consecutive variables Q = O+i = (see Fig.l). In this case the chain 
decomposes into disjunct clusters and Zjv(£) factorizes into independent terms. 

We finally remark that for £ = 1, Z^{Q) in eq.Q can be thought of as a ID version of 
models describing adsorption of hard-molecules, such as, e.g., the so-called "hard-squares" |S] 
or the "hard- hexagons" models 0. 

Periodic CS distribution. - Consider first a situation with a periodic, (with period L), 
CS distribution, such that Q — 6(i,nL + 1) with n = 0, 1, . . . , max(n), where 5(k,m) is the 
Kroneker-delta symbol, 5(k,m) — 1 if k — m and 5{k,m) = otherwise, max(n) = [(N —1)/ L] 
and [x] denotes the integer part of the number x. In this case, the CS density is p = 1/L. 

We have now to distinguish between two situations: namely, when L > 3 and when 
L = 1 or L = 2. In the former case, evidently, the factors (1 — Ci n i n i±i) m eq.Q are non- 
overlapping; Zn(C,) decomposes into elementary three-clusters centered around each catalytic 
site and (possibly) into uncoupled, "free" sites, unaffected by any of the factors (1 — Q n i n i±l)- 
On the other side, for L = 1 or L = 2 we deal with a single cluster spanning the entire chain. 
In fact, the role of L = 1 and L = 2 is, chemically speaking, identical. 

In the case L > 3 the partition function in eq.Q) decomposes into the product 

Z N (C) = zf 3 Z? 2 Zf\ (2) 
where Zx, (K = 1,2,3), are the partition functions for 1-, 2- and 3-clusters, respectively: 
Zx = Yl ^ ll -(l + ^); Z 2 = z ni+n2 {l~n in2 ) - (l+2 2 ); 

{n 1= 0,l} {ni,ri 2 =0,l} 

Z 3 = z ni+n2+n3 (l-n 1 n 2 )(l-n 2 n 3 ) = (l + 3z + z 2 ), (3) 

{rai,ra 2 ,n 3 =0,l} 

while Nk stand for the numbers of such clusters in the iV-chain. Noticing next that A3 = 
max(n) — 5(max(n), (N — 1)/ L), M 2 = 1 + 5(max(n), (N — 1)/L) and using the conservation 
law 3A3 + 2AT 2 + Wi = N, we find that here the pressure per site is given by 

(3P {re9) {L)= lim (\nZ N (Q/N) =p ln(l + 3z + z 2 ) + (1 - 3p) ln(l + z), (L > 3). (4) 

Now, in the periodic case with L = 1 or L = 2 one finds that the partition function of an 
iV-site chain obeys the three-term recursion, Zn — ZjV-1 + zZn-2, whose first three terms 
are given by eqs.Q. Consequently, the Z^ are polynomial functions of the activity z, so that: 

Zn= £ f - +1 )z^z^^F N+2 (l/V-z), (5) 

where Fjss(x) are the Fibonacci polynomials [7] and (^) denote the binomial coefficients. In 
this case, the pressure p( re f) per site for L = 1 or L = 2 turns out to be: 

pp(re g)=ln (Vl + te+ly (L = 1 or L = 2) . (6) 



4 



EUROPHYSICS LETTERS 



By differentiating the results in eqs.lQJ and © with respect to the chemical potential [i, we 
find that in the asymptotic limit z — > oo the mean density n^ re9 ^ of adsorbed particles is 

' 0(\\ L>3, 




- — + 0(^), L = l and L = 2. 

This signifies that for L = 1 or L = 2 and z —> oo the system undergoes an ordering transition. 

Random CS distribution. Annealed disorder. - In this case the disorder-average pressure 
p( ann \p) per site obeys j3P ( ~ ann \p) = limjv^oo (In < Z N (Q > /N) and our aim is to evaluate 
< Zn(() >■ Averaging Zn(Q in eq.Q, we obtain 

<z N (o>= w^-A C(i- P )EL*«y ( 8 ) 

where ^ is the three-site indicator function of the form — (yiirii+i + niUi-i — n,_inirit+i). 
Note now that \Pi always equals zero for unoccupied sites, (rij = 0), and assumes the value — 
1 only for the occupied sites, (n, = 1), which have at least one (or two) occupied neighboring 
sites. Consequently, one has that Yli=i = ■A/+[{ n »}] — ■Ms[{^}] ) where is the 

number of lattice sites on which (in a given realization {n;}) the occupation variable 
assumes the value 1, while -A/i s [{r£i}] is the realization-dependent number of isolated occupied 
sites (elementary cells of the form (0, 1,0)). Hence, < Zjv(C) > m eq.© can be rewritten as 

N N N-N+ + 1 

<Z N (0>= £ (z(l-p)) + (l-p)- m M m (N + \N), (9) 

where M m {N + \N) stands for the number of realizations {n^} that have a fixed N+ and contain 
exactly m elementary cells (0, 1, 0). Using combinatorial arguments, we get [5]: 



Mm {N + \N) = ± u [ N - + 1)l J d ~r( N+ • J [ ) . ,10, 

v +l ; 2m m!(#_ + 1 - m)\ J c r ' 1 - 



itf— JV+ + 1— m 



where C stands for any closed contour which encircles the origin counterclockwise. Substituting 
M m (iV+|iV), ea. l|lU|) . into cq.©, and performing the summations, we are able to determine 
< -Zat(C) >, and hence, also p( arm )(p); in the thermodynamic limit p( an ™)(p) obeys: 



PP {ann \p) = In Uz(l-p)/ 1 - 6v/Qsin (~ arcsin [r/ ' 



(11) 



where 

Next, differentiating eq. ljTJ with respect to /i, one finds the particles' mean density 
n( ann )(p). The resulting expression is rather cumbersome and we present it elsewhere [8]; 
here we merely display the asymptotic behavior of n^ ann \p) in the large- z limit. 

We note that in this limit the forms of rv- ann ' (p) for p < 1 and for p = 1 are quite different, 
which implies that p = 1 is a special point. For p < 1 and z» (1 — p)~ 2 , we obtain 
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while for p = 1 and z — > oo the particles' mean density obeys the expression valid for L = 1 
or L = 2, given in the second line of eq.Q). Note that for p arbitrarily close but not equal to 
unity, n( an ™) = l as z — > oo, while for p = 1 one has n' ann ' = 1/2 as z — > oo. 

This behavior can be understood as follows: In the annealed disorder case, instead of 
averaging lnZjv(C)i we can average Zjv(£) itself. Then, the disorder-averaged pressure is 
defined by the "effective" partition function in eq.JSJ. Here, a strict constraint that no two 
particles can occupy neighboring sites if at least one of them sits on the CS, is replaced by a 
more tolerant condition, which allows for such pairs to be present but a penalty e = 2 ln(l —p) 
is to be paid. For any p < 1, the penalty |e| < oo is finite and one thus expects that for p < 1 
and P/j, ^> |e| the leading behavior is that of the trivial Langmuir adsorption model; hence, 
n (o.nn) ^ z /(l + z) — ► 1 as z — ► oo. On the other hand, for p = 1, e becomes infinitely large 
and can not be compensated by increasing /3fi. 
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Fig. 2 - Particles' mean density n versus the CSs mean density p for different values of z = exp(/3[i). 
Curves with the signs (ann) and (quen) depict the behavior of n in case of annealed and quenched 
random CS distributions, respectively. 



Fig. 3 - Particles' mean density n in the annealed and quenched disorder cases versus the activity z 
for several different values of p. 



Random CS distribution. Quenched disorder. - We finally turn to the most challenging 
situation, namely that of quenched random CS placement. Here we only briefly outline the 
derivation, which is rather lengthy and will be presented in detail elsewhere [S]. 

Consider an TV-chain which contains a fixed number N — N nc of catalytic sites, and hence 
N nc non-catalytic sites, the latter being placed at the positions {X„}, n = 1,2,..., N nc . Next, 
we introduce a set of N nc + 1 intervals connecting consecutive non-catalytic sites, such 
that l n — X n — X n _i (with Xq = 0) and lN ne +i = N + 1 — Xpf ne . That is, the first interval 
reaches from the boundary site i = to the first non-catalytic site, closest to it, the second 
interval goes from this site to the next non-catalytic site and so on, while the interval lN nc +i 
goes from the last non-catalytic site of the chain to the site i = N + 1. 

Now, the logarithm of the partition function in eq.Q, averaged over all realizations of the 
quenched random variable {Ci}, can be formally written as 

JV 

<\nZ N (()>= ^"(l-P^E^^liW). (14) 

JV„ o =0 {Z„} 
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where the inner sum with the sign {l n } extends over all possible solutions of the equation 

h + h + h + --- + lN nc +i=N + l, k>l, (15) 

while ZnHIu}) stands for the partition function corresponding to a given set {l n } of intervals. 
For each set of intervals ZN({l n }) decomposes into smaller clusters, 

Z N ({l n }) = zf l({U|JV) zf 2({i " }|A,) Z? Ml " }lN) . . . Z^ [{l ^ N \ (16) 

where Zk, (K = 1,2,..., N), is the if-cluster partition function, which obeys eq.JSJ) (with N 
replaced by K), while ATr ({ln}\N) denotes the (realization-dependent) number of if-clusters 
in the iV-chain with N nc non-catalytic sites. Taking the logarithm of both sides of ea. (|16fl . 
we obtain 

N 

\nZ N ({l n }) = J2 M K ({ln}\N) \nZ K , (17) 

K=l 

and consequently, < \tiZn(C) > has the form: 

N 

<lnZ N (C)>= E 

wr,n{p)^Z k , (18) 

K=l 

where u>k,n(p) is the statistical weight of the if-clusters in an iV-chain: 

N 

w K<n (p)= E p N ' Nnc 0-~P) Nnc E-MOrJM- (19) 

N nc =o {;„} 

The weights wk,n{p) can be also evaluated combinatorially after some straightforward 
but rather tedious calculations we find that for K ^ 1 and K ^ N: 

wk,n(p) = P {K - 1)j2 (1 - p) iK+3)/2 {2F K (0) + (1 - p)(N -K- 1)F K _ 2 (0) }, (20) 



where <f> = \/p/(l — p) and Fk{4>) are the Fibonacci polynomials |7j, while for K = N one 
has 

w N , N (p) = P N/2 (l-p) N/2 {2F N -i {(j>) +<t>F N _ 2 {<l>)+<t>F N (</.)}. (21) 

Eventually, we arrive at the following exact expression, which determines the disorder-averaged 
pressure p(i uen ) (p) per site in the quenched disorder case (for arbitrary p and z): 

1 1 N 

p P ( q uen) {p) = ^ -(l n Z N (C) ) = lim^ - wk,n (p) In Z K = (1- pf ln(l + z) + 

K=l 

2 „fl _ „\i 



( ,\, / v / T+4i+l\ p(1-p) 2 , , x P(l-p) 
>(5-7p + 3p 2 ) ln(^ J - n 2 P) ln(l + 4z) - _^_^L= x 



V 2 J 2 Vp(4-3p) 

— ; — m+5 



X 

m=0 



2p(l-p) \ / -2p(l-p) 



V^(4 -Zp)- P ) \ vV /p(4-3p)+ P y 
Differentiating ea. ()22jl and turning to the large- z limit, we find that 



lim n(« ne "Vp) = 1 ?—s. (23) 
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This is different from the behavior observed in the annealed disorder case: lim^oo rv- ann ' (p) = 
1 for any p < 1; it also differs from our result for the periodic CS distribution: limz^oo rv- re9 '{L) 
= 1 - p for p < 1/2 and lirn^oo n^ re9 \L = 1 or L = 2) = 1/2 for p = 1/2 and p = 1. Note 
that here, distinct from the annealed disorder case, rS quen ^ (p) does not show any discontinuity 
in the limit z — > oo and p — > 1. The behavior of the particles' mean density in the quenched 
disorder case is also given in Figs. 2 and 3. 

Conclusions. - We have presented here an exact lattice solution describing the equilib- 
rium properties of the heterogeneous catalytically-activated A + A — > reaction in the case 
when the A particles undergo continuous exchanges with a reservoir and react immediately 
if at least one A of a neighboring AA-pair sits on a catalytic site. We have considered three 
possible situations for the CSs placement on the ID lattice: regular, annealed random and 
quenched random. In all three cases we have calculated the disorder-averaged pressure of 
the adsorbate exactly and have presented asymptotic results for the particles' mean density. 
Remarkably, at equal mean densities of the CSs, the asympotic values of the mean densities 
differ in all three cases. We close by noting that the model studied here furnishes another 
example (see, e.g., Refs. [9lll()| ) of a ID Ising-type system with random multi-site interactions 
which admits an exact solution. 
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